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COMPENSATION OF HUMAN VARIABILITY IN PULSE OXIMETRY 
Field of the Invention 

The invention relates generally to. pulse oximeters used to detect blood 
oxygenation. More specifically, the invention relates to a method for taking into 
account human variability, in pulse oximeters. The, invention further relates to; a 
sensor allowing compensation for the inaccuracies caused by human variability, 
the sensor being an integral part of the pulse oximeter. 

\ Background of the Invention 

v Pulse oximetry, is at present the -standard of care for continuous 

'monitoring: of arterial * oxygen saturation (SPO2). Pulse oximeters provider 
instantaneous in-vivoi measurements of -arterial oxygenation, and thereby; ;an 
early warning of arterial: hypoxemia, for example. . 

A pulse oximeter comprises a computerized measuring ' unit ^and a 

» probe attached to the patient, typically to a finger- or ear lobe. The . probe 
includes a light, source for sending an optical signal -through the tissue-and a 

.photo detector for receiving the signal after transmission through the tissue. On-, 
the basis of the. transmitted and received signals, light absorption by the tissue 
can be determined. During- each -cardiac cycle, light absorption by the tissue 

-varies cyclically: 5 : During the diastolic phase,, absorption is caused by venous; 

% blood,, tissue, bone, and pigments, whereas during the systolic phase, there 1 is, 
an increase in absorption, which is caused by the influx of arterial blood into -the'; 
tissue. Pulse oximeters focus the measurement on this arterial, blood portion by 
determining; the difference between the* peak absorption during the- systolic- 
phase and the constant absorption during the diastolic phase. Pulse oximetry is - 

•thus based on the assumption that the pulsatile component of the absorption is: 
due to arterial blood only. 

Light transmission through an ideal absorbing sample is determined by 
the known Lambert-Beer equation as follows: 

(1) 

where= h n is the light intensity entering the sample, lout is the light 
intensity received from the sample, D is the path length through the sample, e is 
the extinction coefficient of the analyte in the sample at a specific wavelength, 
and C is the concentration of the analyte. When l iru D, and e are known, and l 0 m 
is measured, the concentration C can be calculated. 



In pulse oximetry, in order to distinguish between two species of 
hemoglobin, oxyhemoglobin (Hb0 2 ), and deoxyhemoglobin (RHb), absorption 
must be measured at two different wavelengths, i.e. the probe includes two 
different light emitting diodes (LEDs). The wavelength values widely used are^ 
660' nm (red) and 940 nm (infrared), since the said two species of hemoglobin 
have substantially different absorption values at these wavelengths. Each LED 
is illuminated in turn at a frequency which is typically several hundred Hz. 

The accuracy of a pulse oximeter is affected.by several factors. This is 
discussed briefly in the following: 

5 Firstly, the dyshemogtobins which do not ^participate in oxygen 
transport, i.e. methemoglbbin (MetHb) and carboxyhemoglobin (COHb)i absorb, 
light at the wavelengths used; in the measurement. Pulse oximeters are set up 
ito measure oxygen ^saturation . on the assumption that the patient's blood 
composition is ithe -same as r that of a healthy, non-smokinjg • individual. 
Therefore, if these species, of hemoglobin are present in higher concentrations 
- than normal, a pulse oximeter may display erroneous data. 

Secondly, intravenous dyes used : for diagnostic purposes may cause, 
considerable deviation in pulse oximeter readings. However^the effect of these 
dyes is short-lived since the liver purifies blood efficiently. 

Thirdly, coatings like nail polish may in practice impair the accuracy of 
a pulse oximeter, even though the absorption caused by them is constant, not 
pulsatile,, and thus in theory it should not have an effect on the accuracy- 

^Fourthly, the optical signal may be degraded by both noise and motion 
artifacts. One source/ of noise is the ambient light * received by the 
photodetector. Many solutions have been devised ,with the aim of minimizing or 
eliminating the effect of the movement of the patient on the signal,, and the 
ability, of a pulse oximeter to function correctly in the presence of patient motion 
depends on the design of the pulse oximeter. One way of canceling out the 
motion ortofaot *is to use an extra wavelength for this purpose. 

A further factor affecting the accuracy of a pulse oximeter is the 
method used to calibrate the pulse oximeter. Usually the calibration is based on 
extensive empirical studies in which an average calibration curve is determined 
based on a high number of persons. By means of this calibration curve, which 
relates the oxygen saturation of blood to pulse oximeter signals, the average 
difference between the theory and practice (i.e. in-vivo measurements) is taken 
into account. The calibration curve typically maps the measured in-vivo signal 
to a corresponding Sp0 2 value. 



Pulse oximeters, however, can also utilize the Lambert-Beer model for 
calculating the concentrations of the different Hb species. In this method of 
calibration, the measurement signals must first be transformed into signals 
applicable to -the Lambert-Beer model for calculation. This transformation 
5 constitutes the calibration of the pulse oximeter, since it is the step which 
adapts the in-vivo signals to the Lambert-Beer theory, according to which the 
pulse oximeter. is designed to operate. Thus, the calibration curves can also be 
in the form of transformations used to adapt the actual in-vivo measurements to* 
the Lambert-Beer model- 
10 Transformations are discussed for 'example in U.S^Patent 6,104,938, 

which- discloses a calibration method based on the absorption; properties of 
each hemoglobin component, i.e. on the extinction coefficients of- blood r. In this 
-method; the effective extinction coefficients are determined for each light signal 
via a mathematical transformation from the, extinction 'coefficients according to 
}t 15 the LamberfcBeer theory. 

However, each/ patient (i.e. subject of the -measurement) has -a 
^ calibration curve of his or her own, which. deviates from the average calibration, 

- curve calculated on the basis of a high number of patients. This is due to the- 

CH : fact that for each patient the characteristics of the tissue through which light is 

s ^ 20 transmitted deviate from those of an average patient. One drawback of the 

Jjj current pulse -oximeters is -that they are . incapable of taking this human 

t& variability- into account. Human variability here refers to any and all factors, 

causing patient-specific variation in the calibration, curve, including time- 
dependent changes in the calibration curve of a single , patient. As discussed in 
25 the above-mentioned U.S. Patent, subject-dependent variation . can also be 
seen as an effect of a. third substance, such as a third r hemoglobin species in; 
the blood. However, the variation can also be interpreted as a subject- 
dependent change in the calibration curve of the pulse oximeter. 

Without compensation for human variability, the accuracy of current 
30 pulse oximeters is about ±2%Sp02. However, in multi-wavelength applications 
in general, and especially if weak absorbers, such as COHb, are to be 
measured, the human variability represents a much more serious problem. 
Therefore, techniques of compensation for these inaccuracies are called for. 

It is an objective of the invention to bring about a solution by means of 
35 which the effects caused by the tissue of the subject can be taken into account 
when a pulse oximeter is calibrated. In other words, it is an objective of the 
present invention to create a pulse oximeter which can take into account the 
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differences caused by an individual subject as compared to the average 
calibration or transformation curve which the current pulse oximeter relies on. 

A further objective of the invention is to bring about a general-purpose 
solution for the compensation of inaccuracies caused by human variability in 
5 pulse oximetry, a solution which is not limited to the particular general 
calibration method employed in the pulse oximeter, but which can be applied to 
any pulse oximeter regardless of its current built-in calibration method. 

Summary of the Invention 

10 These> and other objectives of 'the invention are accomplished in 

accordance .with the principles' of. 'the * present invention by tprovidingra 
mechanism r-by^-means of which the -.subject-specific deviation in .the tissue- 
induced effects on < the accuracy of the pulse oximeter can be? taken into 
account. 

15 In the; method of the invention, the effect of tissue is ^taken into, 

account and the inaccuracies caused by subject-specific variation in that effect- 
are compensated for. This is implemented by defining a nominal calibration for? 
t-the" -.apparatus- and making off-line: -measurements, also termed "initial^ 
characterization • measurements", in order to define the characteristics which 

20 . describe the conditions under which the nominal calibration has been defined; 
Reference data indicating the characteristics are stored for subsequent on-line* 
measurements- in which -light transmission through the actual tissue is 
measured. (Off-line here refers to measurements performed before the 
apparatus es taken into use, whereas on-line refers to the actual in-vivo 

25 measurements.) Subject-specific -calibration, is then defined based on the 
nominal calibration, the on-line measurements, and the reference data created' 
in connection with the off-line measurements. Thus, the inaccuracies are 
eliminated; by means of on-line measurements, which indicate the effect of the 
tissue. Off-line measurements are used to create the reference data- so that 

30 light transmission measured subsequently through the tissue of a subject can 
be used to correct the nominal calibration for that particular subject. 

Thus in one aspect the invention provides a method for compensating 
for subject-specific variability in an apparatus intended for non-invasively 
determining the amount of at least two light absorbing substances in the blood 

35 of a subject and being provided with emitter means for emitting radiation at a 
minimum of two different wavelengths and with detector means for receiving 
the radiation emitted, the method comprising the steps of 



- calibrating the apparatus using a nominal calibration, 

- carrying out initial characterization measurements, said 
measurements to include the measuring of radiation received by the detector, 

- based on the characterization measurements, establishing nominal 
characteristics describing conditions under which the nominal calibration is 
used, 

- storing, reference data indicating the nominal characteristics 
established, 

- performing in-vivo measurements on a living tissue, wherein radiation 
emitted through the tissue and received by the detector means is measured; 

- based on the in-vivo^ measurements and. the reference data stored, 
determining tissue-induced changes in.the nominal characteristics; and ; 

- :Compensating for subject-specific variation in the^ * in-vivo 
measurements by correcting the nominal calibration on the' basis of the -tissue- 
induced changes. 

In a preferred embodiment of the invention the method : is divided Wn 
two steps so that the first step compensates for the inaccuracies caused ; by 
» tissue-induced- wavelength shift and: the second step compensatesr for the 
inaccuracies caused by interna ^ effects occurring in the tissue. The first step is 
' then* used to -correct the: extinction coefficients -.of- the blood analytes to 1 be 
measured, and the second step isxused to correct the average transformation 
stored in the pulse oximeter 

In a further preferred: embodiment of the invention the effect of the 
temperature is also compensated for in connection with the first step^ 

The method is not^ limited to pulse oximeters explicitly using the 
transformations, but- can . be^ applied • to any pulse oximeter. However, the 
method is preferably appliedfto a pulse oximeter based on a transformation, 
since in a preferred embodiment the method is implemented by carrying out 
changes separately in the transformation and in the extinction coefficients. 

In another aspect, the invention provides an apparatus for non- 
invasive^ determining the amount of at least two light absorbing substances in 
the blood of a subject, the apparatus comprising 

- emitter means for emitting radiation at a minimum of two different 
wavelengths, 

- detector means for receiving said radiation at each of said 
wavelengths and producing at least two electrical output signals, 



- first signal processing means for processing said output signals and 
producing a modulation signal for each wavelength, whereby each modulation 
signal represents the pulsating absorption caused by the arterialized blood of 
the subject, 

- second signal processing means for applying a predetermined 
calibration on said modulation signals, whereby transformed modulation signals 
applicable in the Lambert-Beer model are obtained, 

- memory t means, for storing reference data indicating nominal* 
characteristics under which said predetermined calibration has been applied, 

- first* compensation means, . operatively connected * to the memory 
means, for determining tissue-induced changes in the nominal characteristics, 

- second* compensation means, operatively connected to the first 
compensation-means, for defining a subject-specific calibration by correcting, 
the predeterminedrcalibration on the basis of the tissue-induced changes, and 

- calculation means, responsive to the second compensation ^mearis,. 
for determining said amounts. 

rln a still; further aspect, the invention provides a sensor for collecting ■ 
measurement data for a pulse oximeter intended for non-invasively determining;,, 
the amount of at least two light absorbing; substances in the blood of a subject, 
the sensor comprising- 

- emitter means for emitting radiation at a minimum of two different 
wavelengths, 

- detector means for receiving said: radiation at each of said' 
wavelengths and for. producing at least two electrical output signals, 

- storage; means including reference; data indicating, nominal 
characteristics describing calibration conditions of the^pulse oximeter, said datav 
allowing apparatus connected to the sensor to determine tissue-induced 
changes in the nominal characteristics when radiation is emitted through said 
tissue. 

Preferred embodiments of the invention are discussed in more detail 

below. 

Brief Description of the Drawings 

In the following, the invention and its preferred embodiments are 
described more closely by referring to the examples shown in FIG. 1 to 12 in 
the appended drawings, wherein: 



FIG. 1 illustrates the basic embodiment of a pulse oximeter according 
to the present invention, 

FIG. 2 illustrates the signals utilized in the pulse oximeter of FIG. 1, 

FIG. 3 shows the extinction coefficients of two different species of 
hemoglobin as a function of wavelength, 

> FIG. 4a to 4f illustrate the average transformation curves for two 
different pulse oximeters, 

FIG; 5Js a flow diagram illustrating the prior art calibration method; 

FIG. 6 is a flow^diagram illustrating the general principle according to 
the present invention, 

-FIG. 7 is a flow diagram illustrating compensation for the inaccuracies 
caused by tissue-induced wavelength shift, 

FIG. 8; illustrates »ari- example of ; the transmission : curve of 'human 
tissue, thej curve being employed: in- the compensation -.:of ; ahe inaccuracies 
caused by. tissue-induced wavelength shift, , 

* FIG. 9 illustrates the emitter, characteristics determined for the 

, - • , . 

compensation,- 

>FIG: 10 is a flow diagram illustrating the compensation of the 
- inaccuracies causedby internal effects occurring in the tissue, 

* FIG. 11 illustrates the low frequency baseline. fluctuation! utilized in a 
.preferred -embodiment of the invention; and* 

- FIG. 12 illustrates an embodiment of a sensor according to ? the 
invention. 

Detailed Description of the Invention 

Below; the solution according, to' the invention is discussed with 
reference to a pulse oximeter utilizing the above-mentioned transformations 1 
and; four different wavelengths. As mentioned above, U.S. Patent 6,104,938 
discloses a pulse oximeter utilizing the transformations. 

FIG. 1 is a block diagram of a pulse oximeter utilizing four different 
wavelengths. Light from four different LEDs 10a, 10b, 10c, and 10d; each 
operating at a respective wavelength, passes into patient tissue, such as a 
finger 11. The light propagated through or reflected from the tissue is received 
by a photodetector 12, which converts the optical signal received into an 
electrical signal and feeds it to an input amplifier 13. The amplified signal is 
then supplied to a control unit 14, which carries out calculation of the amount of 
the Hb-derivatives in the blood. The control unit further controls the LED drive 
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15 to alternately activate the LEDs. As mentioned above, each LED is typically 
illuminated several hundred times per second. 

When each LED is illuminated at such a high rate as compared to the 
pulse rate of the patient, the control unit obtains a high number of samples at 
5 each wavelength for each cardiac cycle of the patient. The value of these 
samples (i.e. the amplitude of the received signal) varies according to the 
. cardiac cycle of the patient, the variation being caused by the arterial blood, as 
mentioned above. The^ control unit 14 therefore - utilizes four measurement 
signals, as shown in FIG; 2, each being received at one of the wavelengths. 
1 0 In order for variations in extrinsic factors, such as the brightness of the 

LEDs, sensitivity of the detector, or thickness of-the finger,* to have no effect on 
• the; measurement, each signal 1 : received is normalized by extracting, the AC 
L ^ component oscillating at th£ cardiac rhythm of the patient, and then dividing the 

□ AC component by the DC component of the light transmission or reflection. The 

15 signal thus obtained is independent of the above-mentioned* extrinsic factors. 
Thus in this case the control unit utilizes four normalized signals, which are in 
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M the following denoted with dA i =— — f , where i is the wavelength in question (in \ 
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this basic embodiment of the multi-wavelength, pulse oximeter 1=1,2,3, 4), Ad is- 
p the AC component at wavelength /, and i DC/ is the : DC component at 

fU 20 -wavelength /'. The signals dA are also referred to below as modulation signals. 

J* The modulation signals thus indicate how absorption is affected by the arterial 

p blood of the patient. 

fU The above-described measurement arrangement corresponds to a 

conventional four-wavelength pulse oximeter. The method: of the present- 

25 invention is implemented in the control unit of the pulse oximeter on the basis of 
the four modulation signals described above, i.e. the novelty of the system 
resides within the control unit itself. However, to be able to perform the self- 
calibration in conjunction with each patient, the control unit requires some pre- 
calculated data, which is stored in the memory (M1) of the pulse oximeter. 

30 Instead of being stored in conjunction with the control unit, this data, or at least 
part of it, can also be stored in the sensor part of the pulse oximeter. The 
sensor part, including at least the LEDs and the photo detector, is connected to 
the signal processing part, which includes the control unit. Consequently, 
depending on the overall configuration, the novelty can also reside partly in the 

35 sensor. The operation of the pulse oximeter is discussed in more detail below. 

The theory of pulse oximetry is generally presented as being based on 
the Lambert-Beer Law. According to this theory, light transmission through the 
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tissue at each wavelength is exponentially dependent on the absorbance of the 
tissue (Eq. 1). This theory is generally accepted and established as a good 
model for pulse oximetry. 

Next to be discussed is the theory and formalism on which the method 
5 of the invention is based. 

According to the Lambert-Beer theory and for a system of two 

analytes, the signals described above can be presented as follows: 

dA,=dAx {e? bo > x Hb0 2 + £™ b x RHb) 

dA 2 =dAx (£ f ° 2 x Hb0 2 + £™ b x RHb) 
10 dA 3 =dAx (£? b ° 2 x Hb0 2 + £* Hb 'x RHb) 

dA 4 =dAx(£^- xHb0 2 +£* Hb xRHb) 
RHb=\-Hb0 2 

where dA is a common factor which depends on the absolute values,, 

HbO 

he. inter alia on the total amount of hemoglobin, e, 2 is the .-. extinction 
15 coefficient of oxyhemoglobin at .wavelength i (i=1-4), e^'is the extinctions 
coefficient of deoxyhemoglobin at wavelength i, Hb02 is the concentration-, 
fraction of oxyhemoglobin, and RHb is the concentration fraction of { 
deoxyhemoglobin: 

Using a matrix notation, the above dependencies can be expressed for 
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where dA is the differential change in absorption (i.e. the modulation 
signal) at wavelength Xi, ey is the extinction coefficient of the hemoglobin 
derivative HbXj at wavelength ?u and the constant C accounts for the change of 
25 units to fractional percentages of the concentrations of the analytes HbXj. 

FIG. 3 shows the extinction coefficients (£ Hb ° 2 ar\6 e RHb ) of 

oxyhemoglobin (HbC>2) and deoxyhemoglobin (RHb) as a function of the, 
wavelength. Point P shown in the figure is the isobestic point of oxyhemoglobin 
(Hb02) and deoxyhemoglobin (RHb). The point has the special property that 
30 the modulation signal at the wavelength in question does not depend on the 
respective proportions (relative concentrations) of the hemoglobin species. 
Thus at the wavelength of point P the effect of the relative concentrations of 
oxyhemoglobin and deoxyhemoglobin on the result of the measurement is nil. It 
should be noted, however, that the modulation signal is independent of the 
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relative concentrations only, not of the absolute concentrations. Thus, the 
absolute amount of the hemoglobin species has an effect on the result of the 
measurement. 

As is known, there is a difference between the Lambert-Beer theory 
5 and the practical measurements. The difference is due to the fact that the 
Lambert-Beer theory does not take into account the scattering and non- 
homogeneity of the tissue, whereas the actual extinction coefficients are also 
dependent on the scattering of light caused by the tissue and blood, and on the 
combined effect of absorption and scattering. The larger the proportion of the 
10 attenuation caused by absorption and scattering, the larger is , the correction 
needed^ between the actual and ; the -theoretical (non-scatter) domains. This 
correction .between these two* domains can 'be represented "by the 
-transformation curves' discussed above, by means of which the actual -in-vivo 
□ ■ measurements are mapped to! the Lambert-Beer model. 

1 5 The transformation can be expressed, for example, as follows: 

AtfT* (3), 
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"~ dA 
% where N u =—j- is the modulation ratio (the superscript indicating, the 

tfl dA, 

p domain) in the form of a polynomial function (the subscripts k and / indicating 

/ jm the wavelengths in question), and g is the transformation (the subscrip t- "-1" 
20 denoting the inverse function), 
p ^|l fl Figures 4a to 4f illustrate the average transformation curves measured 

f U for a pulse oximeter, where the two wavelengths for measuring the two species- 

of hemoglobin are 660 nm and 900 nm and the third wavelength is either- 725 
nm or 805: nm. Figures 4a to 4c illustrate the transformation 'curves for a pulse 
25 oximeter with the third wavelength being 725 nm, and Figures 4d to 4f illustrate 
the. transformation curves for a pulse oximeter with the third wavelength being 
805 nm. Each curve shows the Lambert-Beer N k ,i as a function of the in-vivo N M 
at wavelengths k and /. 

FIG. 5 is a flow diagram describing the general measurement principle 
30 described in U.S. Patent 6,104,938. In this method, the above-mentioned 
N%~ vtvo values are first determined from the values measured (step 51). The 

average transformations gw are then used to convert the measured in-vivo 
values to values N^~ B , which can be used in the ideal Lambert-Beer model 

(step 52). Other input values needed for the Lambert-Beer model are also 
35 determined (step 53). In practice these input values are the ideal (nominal) 
extinction coefficients of the analytes to be measured, the extinction coefficients 
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being given for the center wavelengths used in the measurement. The 
converted transformation values and the nominal input values (i.e. nominal 
extinction coefficients); are then used according to the Lambert-Beer model to 
calculate the concentrations of the desired analytes (step 54). Thus in this 
5 approach the in-vivo values measured from the tissue are converted to 

the ideal in-vitro (cuvette) environment, where the ideal oximetry model (i.e. the 
Lambert-Beer model) is applied to yield the desired concentrations. 

In the standard two wavelength pulse oximetry, the prior art technique 
is to map the- modulation ratio W£~ Ww directly to the Sp02 percentage 

10 measured. In this simple case the transformation is not necessary, though the 
transformation technique together with the solution in the Lambert-Beer domain 
can be utilized as well. 

There are two basic ways to determine the average transformation, a 
theoretical approach and an empirical approach. In the empirical approach the 
™ 15 measurements are made in the tissue by taking blood samples and measuring:^ 

\a the- actual" proportions, of the hemoglobin species and then- determining' the' 

\a value of Nu~ B on the basis of the measured proportions. The transformation is 

~ then obtained as the relationship between the values based on the blood; 

5 samples and the values given by empirical measurements as measured by the 

20 pulse oximeter. The theoretical approach, in turn, is:based on a known tissue- 
• model, which takes into account the characteristics of the tissue as referred to 
above, which are ignored* in the Lambert-Beer model. A first value is 
determined for in-vivo N k i by means of the tissue model and a second value on 1 
the basis of the Lambert-Beer model. The tissue parameters of the model are 
25 determined so that the known 2-wavelength , calibration (so-called R-curve) is 
reproduced. Then using these tissue parameters and the wavelength 
dependence of the tissue model, the relation of the in-vivo Nki and the Lambert- 
Beer Nki is extrapolated to other wavelengths in order to obtain the 
transformations at these new wavelengths. Thus in the theoretical approach no 
30 new empirical measurements are made. 

In practice the transformation can be a quadratic equation yielding a 
correction of the order of 20 percent to the measured Nki ,rvv,vo value, for 
example. As discussed below, the transformation data (i.e. the transformation 
curves) are preferably stored in numeric form in the pulse oximeter or the 
35 sensor. The number of transformation curves stored in the pulse oximeter can 
vary, depending on the number of wavelengths used, for example. Typically 
there is a transformation curve for each wavelength pair. 
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As mentioned above, the accuracy of a pulse oximeter utilizing an 
average transformation is not necessarily sufficient, especially if analytes which 
are weak absorbers are to be measured or if two analytes absorb similarly, 
whereby it is difficult to distinguish the said analytes from each other. 
5 In the present invention the accuracy of the pulse oximeter is improved 

by taking into account the subject-specific light transmission through the tissue, 
and changing the values input to the ideal model, i.e. the nominal 
transformation and the nominal extinction coefficients, on the basis of the 
measurement to compensate for the subject-specific changes. 
10 FIG: 6 is a flow diagram illustrating the general principle of the present 

invention. In the method of the invention, initial characterization measurements- 
are first made off-line, preferably at the calibration stage of the pulse oximeter, 
with a nominal wavelength: pulses oximeter sensor (step 61).) Based on»the- 
measurements, nominal characteristics- are established describing, the 
15 conditions under which the pulse oximeter has been calibrated (steps 62): As a * 
result of these steps, reference data are stored (step -63); which describe the- 
l f A calibration conditions of the pulse oximeter. In connection with subsequent in- 

vivo, measurements, the same characteristics are again -estimated and' tissue- 
induced changes in the characteristics are .determined based on the measured 
□ 20 characteristics and the reference data stored (step/64): In in-vivo measurement, 

rU the AT£~ V/TO values are determined from the dA values measured. On the basis 

.of the changes determined, the subject-specific calibration is then determined 
p (step 65) for the in-vivo measurements to be performed by the pulse oximeter 

on the subject. 

25 It is to be noted here that steps 61-63. are off-line steps performed 

either when the pulse oximeter has been calibrated in a known manner or at 
the manufacturing stage of the pulse oximeter sensor when the sensor 
characteristics are determined. After these steps, the nominal transformation 
and the nominal extinction coefficients are known to the pulse oximeter. 

30 Human tissue can influence the accuracy of a pulse oximeter by two 

different mechanisms. First a direct wavelength shift is caused in the LED 
emission due to the filtering effect of the tissue. Namely, on one side of the 
LED center wavelength the Q^porbt i oR may be larger than on the other, 
fjjAfs AfJ^ whereby the center wavelength is effectively shifted towards the region with 
*^y\ 35 smaller absorption. The second mechanism is a subtle one. It arises from the 
fact that the arterial blood is in interaction with the surrounding tissue, which 
can either increase or decrease the effective path length through the arterial 
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blood layer. The first mechanism is in this context termed the external 
mechanism, since it affects factors external to the tissue (wavelength). The 
second mechanism is called the internal mechanism, as it is caused by internal 
factors in the tissue itself: 
5 Therefore, the subject-specific variations in the influence of tissue are 

preferably compensated for by adding two compensating^ processes to the 
above-described prior art mechanism; the first process attends to the subject- 
specific variation in the external mechanism, and the second one to the subject-, 
specific variation in the internal mechanism. The first process preferably 
10 .controls the extinction coefficients to be input to the Lambert-Beer model, while 
the second process preferably, controls the value of the transformation' used to 
transform the modulation ratios N k i" vv ' vo to the Lambert-Beer- model-. N k i L " B . The 
^ linear equations with the unknown analyte concentrations are then solved in the- 

□ Lambert-Beer model, as in. the prior art method. The : degree of -these 

1 5 ■ compensations is determined by DC light transmission through the tissue (the 
vj measured DC signal); measured both in off-line and on-line conditions. . 

M- The first process is illustrated in the flow diagram of FIG^ 7. Initial off- 

line, measurements are again made in order to determine the nominal 
characteristics under which the pulse oximeter has been calibrated (steps 71 , 
P 20 and 72). As a result of these steps,- reference data is created 1 and stored, the 

data indicating the nominal. characteristics (step 73). On the basis of the data 
yj and the in-vivo measurements made^ the tissue-induced changes in the 

O extinction coefficients- are then determined in connection with the in-vivo 

^ measurements (step 74); The new coefficients are then used in the Lambert- 

25 Beer model to solve the concentrations of the analytes (step 75). Thus in this 
embodiment the extinction coefficients are changed; in a subject-specific 
manner in order to improve the accuracy of the pulse oximeter. 

In one preferred embodiment of the invention the compensation for 
any other wavelength shift type interference than tissue-induced interference, 
30 such as a shift caused by high LED temperature, spectral absorption, or 
reflection properties of the sensor cushion or adhesive, is also implemented in 
connection with compensation for the filtering effect of the tissue*. 

In the following the compensations are discussed in more detail by 
disclosing a practical implementation for both compensation steps. The 
35 compensation of subject-variability causing wavelength shift type interference 
(i.e. external mechanism) is discussed first. 
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In the Lambert-Beer model (see Eq. 2) the effective extinctions elj effectIve 
for broadband emitters, such as LEDs, can be calculated as follows: 

ef ective = — JV, (A) * LED. (A(T)) * DET(A) * tissue{X)dX (4), 

where the integration is over the LED emission spectrum LEDi(A(T)), 
5 DET(A) represents the spectral sensitivity of the detector, tissue(x) is the 
spectral transmission of light through the tissue, ss(X) is the spectral extinction of 
the analyte in question, T is the temperature, and W=JLED*DET*tissue*dX 
represents a normalization factor. 

In a preferred, embodiment of the invention, the radiation emitting. 
10 means are: Light Emitting Diodes (LED), but lasers emitting : at one single 
wavelength are also possible: For lasers the effective extinction values are the. 
i extinction-, values at the' laser wavelength, which can depend, however, ontthe- 
D temperature of the emitter component. In the case, of a laser, Eq. 4 is thus hot 

^ needed to calculate the effective extinction value. In the preferred embodiment 

1 5 of the invention the emitter and detector means are located at the tissue site at * 
which the radiation is transmitted through the tissue, but the radiation can also 
^ be conducted 1 to and from the tissue site in a. light conducting fiber or in 

p equivalent conduction- means. In this case Eq. 4 shall also include a term fan? 

fy the spectral' transmission of the radiation conductor. A sensor utilizing light 

20 conducting fibers can be as shown in FIG.* 5 of -the above-mentioned U.S. 
Patent 6,104,938. 

Si** 

fU The extinction coefficients can thus be calculated, according to this 

equation by, determining^ all the above factors, which depend on the actual: 
wavelength values. However, as the task of determining the exact value of Eq. 

25 4 is not possible in connection with a real-time pulse oximeter measurement 1 
using only a few discrete wavelength bands, in practice the result of Eq. 4 has 
to be approximated. The compensation is based on determining nominal 1 
extinction coefficients and approximating their wavelength dependence in 
advance at the factory and using this information in the real measurement 

30 situation to approximate the final subject-specific extinction coefficients. 

The compensation algorithm will now be presented for a 4-wavelength 
pulse oximeter according to FIG. 1, having four LEDs at nominal wavelengths 
of 627nm, 645nm, 670nm, and 870nm. The extinction matrix for RHb (first 
column), Hb02, HbCO, and metHb (last column) and for the above four 

35 wavelengths (627nm on top) is then nominally in L/(mmorcm). 
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^ 1.132 0.1799 0.2734 3.575 ^ 

0.9182 0.1124 0.1337 2.411 

0.7353 0.0885 0.0550 0.5796 

v 0.2071 0.2772 0.010 0.5754, 



(5) 



The above extinction coefficients have been calculated applying Eq. : 4 
at nominal LED drive temperature without the tissue filtering term tissue(X). It 
then represents a nominal extinction matrix for a Sp02 sensor before its 
5 attachment on the tissue site. This extinction matrix is then altered on the basis 
of the measured filtering effect caused by the tissue, when the sensor is 
attached on the site. 

We; next introduce a -parameter called the Functional Light 
Transmission. (FLT)i at a wavelength i, since it is used below in order to make* 
i A 10 al! DCi values (measured at varying. LEDiemission powers at the four discrete* 

wavelengths./) corhparable to each other. Using DC values comparable to each 
other is in practice a prerequisite for unveiling the real effect of the tissue on the. 
measurements/ and the characteristics; of -the tissue. In order to obtain the 
I"* comparable' units, the DC light transmission'for each LED channel (wavelength)* 

15 - is first measured at a certain emitter drive current, and the measured- DC value 
is then reduced in the preamplifier to a detector current, which is normalized to 
Q an emitter current value of 1 mA. When measured without the tissue? in the 

probe, this result is called the* Current Transfer Ratio (CTR) of the probe:?CTR 
\j\ characterizes the sensor design and the efficiency of the light transmission from 

C3 20 the emitters to the detector. It is usually of the order of a few microAmps (of 

^ detector current) per one milliAmp (of LED current). Now the tissue. (e.g«. a 

finger) is inserted into the probe and the CTR is again measured. This result is 
now, called, the Functional Current Transfer Ratio (FCTR) because it is the CTR' 
measured under conditions of the function of the pulse oximeter, i.e. when the 
25 tissue is in place in the probe. The FLTj is then calculated for each emitter 
(wavelength) as follows: 

FLT(emitter#k) = FCTR{\mA - emitter - current) I CTR{\mA - emitter - current) . 

Next the CTR and the FCTR concepts will be linked to the Lambert- 
Beer absorption model and to the actual measured intensities in the pulse. 
30 oximeter. The CTR obviously describes how the external probe design factors, 
such as the color and geometry of the probe, affect the light transmission to the 
detector. On the other hand, the FLT can be associated with the true 
transmission through the tissue in units which are normalized to the emitter 
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efficiency. Therefore, Eq. 1 can be written in a slightly different form, as it is 
often written in transport theory: 

I = I 0 exp(-a * d) = I 0 expC-cr^ * d) exp(-« £t , * d ) (6), 

where d is the tissue; thickness and a is an effective absorption: 
5 coefficient., The above equation has been divided into two components. The 
attenuation factor with omt accounts for all internal absorption effects* such as 
? blood and can be associated to the FLT-value as the FLT equals one when no 
tissue (no internal attenuation) is in the probe, and the factor with oext accounts 
for all external attenuations, such as geometrical factors and multiple surface 
10 reflections without light penetration into the tissue, and can be associated witto 
the CTR of the probe? (The term a" denotes the 'phantom' absorption thickness, 
^parameter for the; external effects.) The term ccext'is mainly., a\Sp02 probe 
design issue which does not influence the measurement accuracy as such, andr 
: thus it need not to be compensated; for by any means. The FLT at wavelength^ 
^ 1 5 can now be defined as: 



□ 



*t FLT * = -j / = expC-a^ * d) = FCTR k /CTRk (7). 
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The FLT thus describes light: attenuation caused by the tissue, and it 
can be related to the DC light transmission in the pulse oximeter. 

It 'is now assumed: that the spectral tissue transmission is as presented 

» 

in 20 in FIG. 8, which shows spectral characteristics of tissue imthe> same units for 

O each wavelength, i.e: FLT as a function of the wavelength; based on an 

empirical measurement, la a continuous real-time Sp02 measurement, the 
transmission is known at 4. distinct wavelength values (the: FLT values derived; 
from the DC values in the pulse oximeter) marked in the figure. At each 
25 wavelength the slope, of the tissue transmission curve can be determined or 
approximated using the four transmission values. The slope then determines 
the change in the tissue transmission in a band of a predetermined width (1 00 
nm in this example) around the center of the LED band. We denote the slopes 
between 627 to 645nm and 645 to 670nm by A and B, respectively. This 
30 definition of the slopes is expressed as: 

FLT(% 2 ) - FLT(A, X ) 
(A 2 ^^ I )*(FZr(A 1 ) + FZr(A 2 ))/2 



and 
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B = 



FLT(A^)-FLT(A 2 ) 



(Z 3 - ^ ) * (FLT(A 2 ) + FLT(Z 3 )) / 2 



400 



where FLT(h) is the measured FLT value determined at wavelength fa 
The estimation of these slopes can be improved by calculating the curvature at 
the center LED (645nm). This curvature (change of the slope/nm) is 

B-A 



curv = 



(A 3 -A ] )/2- 



Finally the expressions are obtained for the slopes s at the three red 
wavelengths using^ and B as parameters: 

(A + B)/2 
B + curv*^-^)/! 



(8). 



where the slope at the IR wavelength has been estimated tto be 
constant as it cannot be determined by the other LEDs. If we had had another 
LED; at about 800-1 OOOrim range, for example, it could have been used for, the - 
estimation of the IR slope. Because the extinction curves are very flat at 870nm- 
and' the ; transmission is usually rather high, the; tissue prefilter cannot alter the 
effective extinction^ coefficient from its nominal- value significantly. The 
approximation of a constant transmission slope is thus considered sufficient. . 

In principle these slopes could be inserted in Eq. 4 in order to integrate, 
the new true, values for the extinction coefficients. However, this is impractical ' 
to do in real-time; so. a simpler algorithm is presented below. We first calculate 
off-line the relative changes of extinction coefficients for each analyte of the 
system using: Eq. 4 and assuming that the value of the slope equals; a 
predetermined' value, which is one in this example. This calculation results in a 
shift matrix of Eq. 9: 

f 0.975 0.942 0.940 0.999^ 
0.984 0.966 0.920 0.916 
0.963 0.986 0.896 0.789 
k 1.01 1.03 0.903 1.05 , 



Tissue SLOPE ~ x = 
l issue SHIFT — 



(9), 



where the effective extinction of Hb02 at 645nm is 0.966 times the 
original value, and the effective extinction of HbCO at 670nm is 0.896 times the 
original value, for example. Thus, the matrix of Eq. 9 defines the relative 
changes caused by the tissue, assuming that the slope of the tissue 
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transmission curve equals one. During in-vivo measurement, the slope is 
estimated using the DC values. The ratio of the slopes then indicates the 
relative change of a coefficient. In other words, if the relative change 
determined off-line for a slope value of X1 is r, the relative change for a 
measured slope value of X2 is then rx(X2/X1). The relative changes are 
different for different analytes since the extinction coefficients of the different 
analytes behave differently as a function of wavelength. 

In real-time the effective extinction coefficients can thus be calculated 
as follows: 

K Eff =E* <8>(1 + S •(Tissue™™ -1)) .(10), 

where S denotes the column array in Eq. 8 and the matrix 
multiplications are performed element by element (<g>) or element by row (•), 
respectively. 

If the LEDs are not driven at the nominal drive currents, their effective 
wavelength may also be' shifted by : the temperature change at- the. LED*p-n 
junction. The wavelength shift . induced by temperature is typically about 0.1-0.2 
nm/°C which is significant if the drive currents are high, as is usually the case at 
wavelengths shorter than 660 hm. Thus, the extinction matrix of Eq. 10 must 
also be compensated for in varying LED drive conditions. 

There are many ways to find out the temperature of the LED :p-n 
junction. One alternative is to add a temperature sensor on the LED/ substrate 
and use the reading of the sensor for the compensation of all LED emission 
wavelengths. Though the junction temperature follows the substrate 
temperature according to : some empirical heat conduction model, the method 
may be unreliable because the LED chip contact to the substrate and the 
internal heat conductivity may vary considerably. A better way is therefore to 
determine the junction temperature directly from the forward voltage drop of the 
LEDf junction. The junction has typical diode characteristics, which can be 
determined off-line for each LED separately after assembling the LEDs on the 
substrate. It is even possible to measure, with an optical spectrometer the shift 
of the emission as a function of the LED forward voltage. Relating the 
wavelength shift to the forward voltage assumes that the forward voltage is 
measured during the operation of the pulse oximeter. The circuit board of the 
pulse oximeter should thus preferably have means for performing the forward 
voltage measurement. But if it does not, the LED emission shift can be 
calibrated against the temperature sensor at the substrate. The LED 
manufacturer specifications for the temperature shift can then be used to 
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calculate the corresponding wavelength shift. Still another practical 
compensation for the emitter temperature changes is to map empirically the 
relationship of the emitter drive current to the observed wavelength shift and to 
use this information to adjust the on-line extinction coefficients for the sensor. 

A method for temperature compensation of the LED emission is now 
presented, assuming that-the LED forward voltage is measured on the circuit' 
board. The wavelength shifts can then be calculated as follows 



(11), 



where the shift coefficients are values determined empirically -in 
advance and aVi are the' measured changes of the forward voltage drops.r.For- 
the 627-645-670-870nm > LEDs of; the sensor, the k-values are 0.06 nm/mV, 

•'0.06rim/mV, 0.09nm/mV„and 011 nm/mv? respectively. 

As in the compensation discussed above relating to tissue filtering, it is 

-practical to first calculate the change imthe extinction coefficients off-line* for a 
certain fixed* wavelength shift; In this example the relative changes *of the> 
extinction coefficients are calculated, as in Eq. 9, for a 5 nm wavelength shift for 
each of the four hemoglobin derivatives. The following shift matrix is ■ then 
obtained: 

T0.919 0.820 0.798 0.974^ 
0.963 0.926 0.82> 0.794 
0.941 0.983 0.855 0.725 
1.0 1.01 0.963 1.02 



Tem P%HiFr m 



(12). 



During in-vivo measurement the relative changes are then calculated 
based on the measured wavelength shift. The ratio of the wavelength shifts 
then indicates the relative change of a coefficient caused by temperature. In: 
other words, if the relative change calculated for a wavelength shift of Y1 is r, 
the relative change for the measured (on-line) wavelength shift of Y2 is 
rx(Y2/Y1). The relative changes are different for different analytes, since the 
extinction coefficients of the different analytes behave differently as a function 
of wavelength. 

The temperature compensated extinction coefficients are thus: 
E^p =E° <8>(l + (AA/5/im).(remp^ -1)) (13), 
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where AX is the array in Eq. 11. As mentioned earlier, AX can also be 
estimated by reading the temperature indicated by the temperature sensor on 
the LED substrate or by measuring the LED drive current and using the 
mapping of the current to the wavelength shift. 
5 The compensation of the variability causing wavelength shift type 

interference can now be summed up as follows: 

E Eff =E° ®(1 + S •(Tissue™™' 1 -\))0(1 + (AX/ 5nm)»(Temp^ m -1)) (14). 

As mentioned above, the practical compensations can be divided into 
steps performed off-line prior to the actual measurements, for example in the 
10 factory at the manufacturing stage, of the pulse oximeter sensor, and to steps 
.performed on-line, i.e. in connection with the actual in-vivo measurements. 
Using-the above practical example, the off-line j compensation steps are^s 
follows. 

1. The; spectral characteristics of the emitter/detector system^ is 
15 established.Mn other words, the LEDs are characterized for their light emission 
(the emission, as a function of wavelength) and; the detector for its spectral 
sensitivity. This step thus includes determination of the characteristics of the 
curve shown in FIG. 9, i.e. the received intensity as a function of wavelength (at 
least around, the wavelengths used). The light transmission from the light 
20 emitter to the light detector is measured without human tissue, i.e. the CTR is 
determined - in^ the;; sensor in a fixed setup: mimicking the actual use* of ithe 
sensor. For clip-type sensors this , is usually, the Probe Off position of. the 
sensor. The step also includes determination of the center wavelength of each 
LED. 

25 2. The effective extinction coefficients are determined without the 

tissue term. Thus in this step Eq. 4 is used without the tissue term (tissue(x)) to 
form a matrix Eki° according to Eq. 5. 

3. The relative changes in the effective extinction coefficients due to 
tissue filter effect are determined: Here Eq. 4 is used assuming that the slope of 

30 tissue(x) equals a predetermined fixed value at each wavelength. In other 
words, the shift matrix of Eq. 9 is determined. 

4. The relative changes in the effective extinction coefficients due to 
wavelength shift caused by changes in temperature are determined. In other 
words, the matrix of Eq. 12 is determined, which indicates the relative changes 

35 for a wavelength shift of a predetermined value. 

5. If the LED forward voltage method is used, the LED forward 
voltages are characterized at a typical drive current for small ambient 
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temperature changes. The temperature coefficients ki in Eq. 11 and the 
nominal forward voltage drops at nominal temperature are determined for each 
LED. 

6. The information obtained is saved in a memory unit in the sensor or 
5 in the control unit, or the corresponding information is otherwise made available 
to the pulse, oximeter, for example, by using codes, such as sensor 
identification numbers, which indicate the values of the information. 

As is obvious from the above, steps 1 to 6 include performing initial 
measurements for the compensation, said measurements to include measuring 
10 the light transmission of the apparatus, establishing: nominal DC transmission 
characteristics of the ■ apparatus on the basis of the measurements, and for 
subsequent- in-vivo measurements storing reference data which indicate the 
iA transmission characteristics established. 

□ Thus in the off-line measurements, the value of Eq. 4 is determined 
£3 15 using nominal values, and the nominal extinction matrix, is formed*. The 

Si 

vj apparatus -is also provided with the; data needed in the subsequent- on-line 

^ compensation steps for calculating the changes in the factors included in EqT4. 

^ In the on-line steps: the said changes are determined and: a new extinction- 

matrix is formed, whereby the new extinction values are such that the external 
p 20 .. effects are compensated for. 

To sum up, after the off-line steps the pulse oximeter stores the 
matrices according^ to equations 5, 9, and 12 and the values of the shift 

□ coefficients ki . In addition to this; the oximeter stores the CTR values and the 
I'M center wavelengths corresponding to these values. 

25 Following the above practical algorithm, the on-line compensation 

steps utilizing the information stored during the off-line steps are as follows: 

1 . The forward voltage drops (aV) are determined for each LED. In this 
step, the change in the forward voltage relative to the sensor nominal values 
are calculated, the nominal values having been stored in the sensor memory 

30 unit or in the control unit memory. If the pulse oximeter does not have forward 
voltage measuring means, the temperature of each LED is estimated by 
reading the temperature indicated by a sensor on the LED substrate, and either 
the manufacturer specifications or empirical data for corresponding wavelength 
changes or a look-up table mapping the emitter drive current to the center 

35 wavelength shift is used. 

2. The effective extinction matrix for the temperature compensation is 
determined. This is calculated as illustrated in Eq. 11-13, i.e. using the 
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wavelength shifts determined according to Eq. 11 and calculating the matrix of 
Eq. 13 using the matrices of Eq. 5 and 12 stored earlier. 

3. The DC light transmission for each LED channel (wavelength) is 
measured, and the value measured is normalized to an emitter current value of 
1 mA. The result is the FCTR of the sensor. An estimate for the FLT is then 
calculated for each emitter (wavelength). In this connection, the FLT values are 
for calculating the slopes (A and B). In other words, all DC values are 
normalizedMn relation to the 1 mA emitter current in order to make all values 
comparable to one another. 

5. The effective Lambert-Beer extinction coefficients are determined 
using Eq. 14 

6. The extinction coefficients obtained in step 5 are used for solving 
the analyte concentrations in the Lambert-Beer model. 

5 During: the above on-line steps \-6 in-vivo measurements, are- 
performed, wherein* the- DC component of the; radiation emitted through? the 
tissue and received by the detector is measured; tissue-induced changes in the 
transmission characteristics are determined based on the in-vivo DC 
component and the transmission characteristics stored, and on the basis of the 
tissue-induced; changes the subject-specific variation in the in-vivo 
measurement is compensated for. 

The above off-line and on-line steps compensate for the non-ideal 
characteristics, of the broadband emitters or for external effects on the light 
source emission spectra. However, they do not compensate for the variation in 
the absorption and- scattering interplay in the tissue; i.e. the internal effects, 
which equally influence a single line laser emitter and a broadband* LED 
emitter. Lasers also show shifts in the emission line wavelength as a function of: 
the temperature. Therefore, the lasers are compensated for the temperature 
and internal tissue, effects, but not for the pre-filter tissue-induced spectral 
shifts. 

The compensation step attending to the variations in the internal 
effects is illustrated in the flow diagram of FIG. 10. Initial measurements are 
again made in order to determine the nominal DC transmission characteristics 
of the pulse oximeter (steps 101 and 102). As a result of these steps, reference 
data are created and stored, and the data indicating the transmission 
characteristics are determined (step 103). On the basis of the data and the in- 
vivo measurements made, the tissue-induced change in the average 
transformation is then determined in connection with the in-vivo measurements 
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(step 104). The new transformation is determined, and finally the 
concentrations of the analytes are solved in the Lambert-Beer model (step 
105). Thus in this embodiment the average transformation stored in the pulse 
oximeter is changed in a subject-specific manner in order to improve the 
5 accuracy of the pulse oximeter. As shown below, in this case during the off-line 
steps 101-103 partly different data is to be created and stored as compared to. 
the first compensation step, since in this case the DC light transmission 
characteristics are used to create a subject-specific transformation. 

A practical implementation of the second compensating step is now 
10 discussed by introducing a new variable called "path length multiplier", since 
this will provide an easy way of understanding the technique in accordance with 
the invention. 

As mentioned above, the purpose of the invention is to improverthe 
t~ accuracy of a pulse oximeter in situations in which the blood volume, the red 

□ 15 bloods cel^ density or the hematocrit;- the total hemoglobin (g/dl), the division 

between- the arterial and venous blood compartment volumes, and the arterial- 
venous saturation difference vary and, produce human variability, which 
worsens the accuracy of the SpQ2 measurement. It is also the purpose' of the 
£H invention to compensate for the effect of skin pigmentation (dark skin), which in 

p 20 part can. 'be considered to belong to the tissue prefilter category of 

ry compensations, 1 but which also influences via modifying the path length 

multiplier. This modification is especially important for Sp02 ear sensors, which 
p are attached to a very thin tissue part (of about the same thickness 

ry corresponding to the diffusion constant in human tissue). 

25 The interdependence of the above-described transformation and the 

path length multiplier is first illustrated by considering the photon path lengths* 
through a single layer of artery blood and examining how the scattering and 
absorption affect it. It is postulated here that multiple scattering effectively 
increases the photon path length through the artery and that the absorption of 
30 the surrounding tissue effectively decreases it. In this way the artery and tissue 
are in interaction with each other. To derive a mathematical formulation of this 
relationship, the known Kubelka-Munk two-flux model can be used. This model 
defines an absorption probability K as follows: 
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where Za is the macroscopic absorption cross-section of the media and 
dl is the true average photon path length through the scattering and absorbing 
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medium of infinitesimal layer thickness dz. The term <dl/dz>=K/£a is a path 
length multiplier (plm) which enhances the arterial blood absorption from that of 
the Lambert-Beer non-scatter value because of the multiple scattering in the 
surrounding medium. 

The idea of the path length multiplier is applied to the Lambert-Beer 
formulation of 2-x pulse oximetry. The ratio of the change in absorption at the 
two probe wavelengths is defined as: 

~dA k 



in— vivo 
kl 



(16). 



where ^ is the arterial (non-scatter) absorption coefficient at 
.wavelength / and d, is the effective true optical path length. The transformation 
is defined by substituting equation 15 with dl = d, in equation ,16: 



N. 



in-vivo 



•dz 



J * 



k_ + *r ideal 
Iy kl 



(17), 



where the ideal Lambert-Beer model is used for Nk! deal = ^ta/^J , and 
where the layer thickness dz is the same for all wavelengths (k, /). Equation 17 
now represents the transformation (g M f from Nki in_vivo f i.e; from the measured 
value, to Nki L " B , which is the ratio of differential absorptions, that would- be 
measured if the measurement system were the ideal cuvette system of the 
Lambert-Beer model. For the transformation the equation below is obtained: 



Ski = 




(18). 



Thus the transforming quantity is a ratio of path length multipliers 
measured at two different wavelengths (fc and /). The dependence of the 
function gw thus refers to the absorption density of the scattering tissue in the 
surrounding of the infinitesimal arterial layer dz including the layer itself. This 
essentially means that the transformation does not require knowledge of the 
analyte composition in the arterial blood, but refers rather to the macroscopic 
light absorption, i.e. transmission through the tissue part under the sensor. That 
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is in the language of pulse oximetry the DC component of the light 
transmission. This is utilized in the compensation of the invention. 

Modifying Eq.1 and leaving the attenuation of the probe design 
factors (i.e. CTR values) out of consideration, the relationship of the DC light 
transmission through the tissue and the path length multiplier can be presented 
as follows: 

1^=1^=1^^ (19), 

where D is the actual path length through the sample, D1 is the 
shortest path length through the sample (i.e. the thickness of the sample), and £ 
is* the ideal extinction coefficient of the analyte\ Here the W//n . can . be 
associated with the FLT at the wavelength in question. P/m thus describes the- 
internal attenuation factors in the tissue and, in, particular, the enhancement of 
the absorbancy relative to the ideal cuvette absorption. 

In nominal conditions, the path length multiplier has a certain nominal, 
value pirn 0 (where the superscript '(X refers to. the nominal value). This nominal , 
value can be determined in the factory at the manufacturing stage of the* pulse/, 
oximeter. When the DC component is measured again in connection- with- in- 
vivo measurement, the change in the p/m from the nominal value can be used 
to determine the change in the average transformation. 

The term otint in Eq. 6 can be expressed with the help of the path length 
multiplier in the Lambert-Beer model as 

where 2a accounts for all internal absorption sources and is defined in 
the. non-scatter Lambert-Beer domain. The FLT at wavelength k can then be 
written as follows: 

FLT k = exp(-a int * d) = exp(-plm *£ a *d) (20). 

We then ratio the logarithms of the FLTs at two wavelengths k and /, 
which results in Eq. 21: 

where g M =plm k /plmi is the transformation between the Lambert-Beer 
and in-vivo modulation ratios according to Eq. 18, and in which the internal 
absorbing tissue compartments are venous and arterial blood with volume 
fractions f v and f a and with absorption coefficients ^ and ^ determined in the 
Lambert-Beer domain, respectively. In the last expression we have used for the 
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venous volume fraction the relationship f v = 1-f a - As the arterial volume fraction 
is always smaller than the venous volume fraction and as the arterial-venous 
absorption difference is always smaller than the venous absorption, the 
dominating factor in the last term is jjf/jjJ , i.e. the venous saturation Sv02. 

Thus the changes in the FLT and Sv02 from their nominal values 
provide the compensation needed for estimating the correct transformation 
function . g*/. We can then finally write for the relative change of the 
transformation function g«: 

\og(FLT k ) 



\of>(FLT,) 

\o%{FLT k )° 

g^ = / \og(FLT t )° 

gl F^SvO^SaO^fW 

/F^SvO^SaOlJS 



' (22), 



.where the function Fm represents the ratio term " , a . 7- 7- in 

;Eq. 21. and the superscript 0 represents the nominal values of the nominal 
calibration function gw° , which is on average true for a large patient population. 
In fact, the log(FLT) and the Fkt compensation terms account for quite different 
human variability factors in the tissue: whereas Fkt mainly tracks the changes of 
the arterial venous saturation difference, in particular Sv02, the log(FLT)" 
reflects the changes} in the total absorption of the tissue, i.e. in the total blood 
volume and the total hemoglobin or hematocrit, which are not seen in Fkt at all. 
In practice, the largest corrections to the transformation function are due to the 
log(FLT) and Fkt is less important. 

The FLT in Eq. 22 is easily obtained at the two wavelengths k and /, as 
has been described earlier in Eq. 7. The function Fkt represents the ratio of the 
absorption coefficients (in the Lambert-Beer non-scatter model) of the whole 
tissue at these same two wavelengths, i.e. it represents the internal color of the 
tissue. This internal absorption ratio can be measured by examining the low 
frequency baseline fluctuations of the plethysmographic wave signal. These 
fluctuations are caused by the low frequency changes (usually of respiration 
origin) in the blood volume or in the blood volume distribution of the tissue. 
Similarly, since the arterial color (= R-ratio) is defined as the ratio of the arterial 
absorption coefficients, function Fkt can be calculated as: 

~ _ /«(/*« ->"v) + _ .-l/ ^^C^ _ -1 s xr baseline \ ,™ 
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where AC is the amplitude (or the instantaneous slope) of the low 
frequency baseline fluctuation, instead of the heart pulse amplitude of the 
plethysmographic wave, and DC is the DC light transmission at that particular 
wavelength. Because, the effective tissue color is mainly determined by the 
5 venous blood, function Fid can be approximated as the arterial modulation ratio 
calculated for the venous saturation, which is usually about SaO2-10% i.e. 
Fki=R(SvO2=SaO2-10%). 

If the venous saturation is determined by venous blood samples and 
the arterial saturation by the arterial blood samples, the function Fw can be 
10 calculated using-- the real blood values (with the assumption that the 
corresponding blood compartment volumes are f a =0.25 and fv=0.75). 

The second compensation step for human variability .can then ' be 
^ performed in the following way: 

Q 1 . The nominal light transmission through a finger or ear lobe or its. 

P 15 approximation at the distinct nominal wavelengths of the sensor, is determined* 

%i for a population of patients or volunteers or even for only one single volunteer, 

\* on whom the nominal calibration (the original standard pulse oximeter 

, calibration or the; actual multi, wavelength measurements using this calibration 

\o % {FLT k )° 



5_ directly or indirectly) was performed; this gives ^ 0 ^ FLT y as a function of N w 

20 ^ at these wavelengths (/c, /), ' e as a function of the correct Sp02 in practice. 



in 



Thus, the FLT° values are obtained for each wavelength. This information is 
o stored in the sensor memory unit (or in the control unit). 

I" 2. Function F° can be determined in two alternative ways: 

2.1 In the above measurement using the baseline ^fluctuations, 
25 A/ w tese/ ' ne is calculated . and transformed to the Lambert-Beer model by the 
nominal transformation (g°ki)' 1 for the population of patients or volunteers or 
even for only one single volunteer, on whom the nominal calibration (the 
original standard pulse oximeter calibration or the actual multi-wavelength 
measurements using this calibration directly or indirectly) was performed; this is 
30 tabulated as a function of A/w n " wv0 at these wavelengths (k, /). These data are 
stored in the sensor memory unit (or in the control unit). 

2.2. In the above measurement venous and arterial blood samples are 
taken from a position close to the sensor site and analyzed for RHb, HDO2, 
HbCO, and metHb. The absorption coefficients ^ and ^ are then calculated 
35 using the measured analyte fractions. The arterial volume fraction f a is then 
estimated. Usually it is sufficient to approximate that f a is equal to 0.25. The 
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functions Fw are calculated in the Lambert-Beer Model using the venous and 
arterial volume fractions and 

MatV = RHb * MoT + HbOl * v"t° 2 * HbCO * Ma,f° + metHb * /C"* . 

3. As a result of steps 1 and 2, we have an off-line determined look-up 
table for each wavelength pair, in which the following nominal information is 
tabulated: 



k 


I 


Nk ,n-v,vo 


„ 0 

gki 


Log(FLTk)/log(FLTi) 


_ -1 /KI baseline* 

gw x(N k i ) 


Fki 

















Table!. 

and where also the nominal transformation g° is presented . and 
Nki basel,ne = gki°x(Fki) .'Only one of the two last columns is necessary, depending;: 
on the way the values of function Fare determined. It is also to be noted that 
the ratio is not necessary in column 5, but that it is enough to store the; FLT 
values from which the ratio of their logarithms can be calculated. 

As is obvious from the above, steps 1-3 again include measuring the 
. DC light transmission of the apparatus, establishing nominal DC transmission 
characteristics for. the apparatus on the basis of the measurements, and storing; 
reference data for subsequent in-vivo measurements, the data indicating the; 
transmission characteristics established: 

After the above off-line steps, which can be performed at the 
manufacturing stage of the pulse oximeter, the following steps are performed in 
connection with the on-line (i.e. in-vivo) measurement. 

4. The actual spectral shape of the light transmission through a finger 

or ear lobe at the distinct wavelengths of the sensor is determined, and 

\og(FLT k ) 



\og(FLT t ) 
stage. 



is calculated for each wavelength pair similarly as in the off-line 



5. In Table 1 in the row where the measured modulation ratio at the 



wavelengths k and / equals N k " 



the nominal 



\og(FLT k ) c 
\og{FLT t y 



is read. The 



togCFXr.) /\og(FLT.)° 
correction factor - — ,^ r ^, / - — is then calculated. 
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6. In Table 1 in the row, where the measured modulation ratio at the 
wavelengths k and / equals Nki 1 ""™ the nominal F*/° is read from either of the 
last two columns. 

7. The correction factor F*/Fw° is determined in one of the following 
5 two ways: 

7.1 Column g' 1 x (h/ ,aseSne ): By using the baseline fluctuations of the 
measured plethysmographic wave and by using Eq. 23 with g'\ the function Fw 
is determined. The correction factor F*/Fw° is then calculated. For determining 
A/ tese/ ' ne and its changes, the amplitudes of the signal can be used, as is 

1 0 normally done for a modulation ratio N. 

7.2 Column F«: The blood analytes RHb, Hb0 2 , HbCO, and metHb are 
solved in the Lambert-Beer model using the nominal transformation g°: This is 
the first approximation for the analytes. The absorption coefficients and; ^ v 

□ are calculated using; the measuredv analyte fractions in the arterial absorption 

P 15 and approximating the absorption coefficient in the venous blood by using the 

measured -dyshemoglobin fractions and setting Hb02 yena =Hb02-1 0% and 
J* RHb vena =RHb+10%, where f a = 0.25 is assumed. The functions F k i are 

calculated in the Lambert-Beer Model using the equation: 

MatV = RHb * M ™ b + Hb02 * ju"t° 2 * HbCO * m™ C ° + rnetHb * . 
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fU 20 A more accurate estimate of F M can be obtained by iteration of new 

J* analyte fractions for the new corrected- transformation after step 8. 

q 8. A new transformation g M is calculated using Eq- 22, and the new. 

fU transformation is used for solving the analyte concentrations in the Lambert- 

Beer model. 

25 The pre-calculated data utilized by the pulse oximeter can be stored in 

the sensor part of the pulse oximeter, whereby the same sensor can be 
attached to different pulse oximeter housings. FIG. 12 illustrates the general 
structure of a sensor according to the invention, the detailed configuration of 
the sensor being dependent on which information is stored in the sensor and 

30 which in the signal processing part; and also on the amount of the calculation 
appropriate in the signal processing part. 

Nevertheless, a sensor according to the invention includes the light 
sources (10a-10c) and the photo detector (12), the light sources being adapted 
to emit at two or more wavelengths. In addition, the sensor includes a data 

35 storage unit M2 for storing the data on the basis of which the signal processing 
part can perform the above-described calibration. The information necessary 
for the above compensations is shown in the figure. For the first compensation 
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process the pulse oximeter needs the k-values, the above-mentioned three 
matrices, i.e. the nominal extinction matrix (Eq. 5), the shift matrices (Eqs. 9 
and 12), and the CTR/wavelength pairs. For the second compensation process, 
in turn, the pulse oximeter needs the information stored- in Table 1 and the 
CTR/wavelength pairs. As mentioned above, at least part of this data 
determined prior to the use of the device for in-vivo measurements can also be 
stored in the: control unit part of the pulse oximeter. The apparatus further 
preferably, includes means 120 for measuring the forward voltage of the p-n- 
junction of each LED; as discussed above. 

As can be seen from the above, the method of the invention is based 
on the DC transmission of light. By means of the off-line DC measurements, 
reference data is first Created. During subsequent in-vivo DC measurements, 
*the reference data is then utilized to filter out human variability, from : the in-vivo 
measurement. 

Although the method in accordance, with the invention has-been 
discussed in ; connection with a four wavelength pulse oximeter, it can/also be 
employed in a basic; two wavelength. pulse oximeter. However, the method is 
more beneficial in a multi-wavelength pulse oximeter where the number of 
analytes to be measured is greater than two. 

In the case of a two-wavelength pulse oximeter, the simplest way to 
apply, the compensation is first to formulate -the calibration of the two- 
wavelength oximeter as a first step using only one transformation function g' 1 
(e.g. at wavelengths 660nm and 900 nm)'and a second step using a two-times- 
two extinction matrix s for these wavelengths and for the two analytes RHb an& 
Hb02. The compensation procedures are then identical to the ones presented 1 
in the above multi-wavelength method. If the calibration of the two-wavelength^ 
pulse oximeter is done in the normal way using a direct mapping of the in-vivo 
measured R-ratio (= Neeo-goo) to the Sp02 percentage, the compensation steps; 
could for example, be as follows: The wavelength shifts from the nominal LED/ 
center wavelength values to a change in the Sp02 value can first be coded. 
The wavelength shifts are determined for the temperature component as 
described in the above multi-wavelength method and for the tissue component 
by mapping at the two wavelengths the change in the FLT ratio from its nominal 
value in the calibration conditions to a change in the Sp02 value from the 
nominal calibration Sp02. The tissue wavelength shift cannot be estimated as 
accurately as in the multi-wavelength oximeter, but sufficient compensation to 
the tissue prefilter can still be obtained and the accuracy of the pulse oximeter 
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can be improved. The last compensation step also includes the compensation 
for the internal tissue variability, which is summed with the prefi Iter effect. 

Although the invention has been described above with reference to the 
examples shown in the appended: drawings, it is obvious that the invention is : 
not limited to these, but may be modified by those skilled in the art without 
departing from the scope and spirit of the invention. For example; instead of 
transformation; any other quantity by which the pulse oximeter can correct the 
average calibration known to it can be used to eliminate human variability. 

The invention has also been described with reference to pulse 
oximeters for analytes which are in the blood of a subject. The invention, 
however, can also be applied at different wavelength ranges, e.g: at around 1.5 
for glucose, at which similar compensation means are called for. The , other 
substances in the tissue modify the effective extinction of the glucose; because . 
• they alter the path length multiplier at this wavelength. Similarly, the tissue 
prefilter and temperature effects are taken into account. 



